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ABSTRACT 

Particle heating in advection-dominated accretion flows (ADAFs) by nonlinear MHD (Alfvenic) tur- 
bulence is investigated. Such turbulence with highly- fluctuating magnetic fields, B ~ Bq, is believed to 
be naturally produced by the magnetic shearing instability near the nonlinear saturation. It is shown 
that the energy is dissipated in the parallel cascade, which occurs due to nonlinear compressibility of 
high-amplitude turbulence, and predominantly heats protons, but not electrons. The conservative limit 
on the electron-to-proton heating fraction is 6 <J few x 1CU 2 . 

Subject headings: accretion, accretion discs — MHD — plasmas — turbulence — waves 



1. INTRODUCTION 



The Advection-Dominated Accretion Flows (ADAFs) 
are a class of hot, optically th in accretion solutions 
(llchimaru 19771; iRccs et al. 1982t jNarayan fc Yi 1995 



Abramowitz et al. 1995 | ; see e.g., Narayan. Mahadevanj 



& Quataert 1998 for review) which describe quite well 



the spectral characteristics of a number of low luminos- 
ity a ccreting black hole systems, e.g. , black hole bina- 



Narayan, Yi, fc Mahadcvan 19951; iManmoto et al. 1997 



rics (Narayan, Barret, fc McClintock 1997 ; Eamcury et 
al. 19971) and low luminosity galac t ic nuclei (e.g., Sgr A 



Mahadcvan 19981; [Narayan et al. 1998|; NGC 425 8: |Lasot 
et al. 1996 1 ; I Cammic, Narayan, fc Blandford 1999 1 ; M87 and 
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model with wind outflows ( Blandford fc Bcgelman 1999| ). 
Because of the lower accretion rates (due to the mass loss) 
at the inner regions, where most of the radiation is pro- 
duced, the density of the inflowing gas and the magnetic 
field strength decrease (assuming constant (3) and result in 
a much low er radiative efficiency of the flow. It was demon- 
strated by Quataert fc Narayan (1999| ) that the values of 
S as high as S ~ 0.3 are still plausible in the advection- 
dominated flows with strong winds. The question of turbu- 
lent particle heating and non-thermal coupling of electrons 
and protons i n hot accretion flows has been addressed in 
several works (Quataert 1998 ; Quataert fc Gruzinov 1999 



i i ^ — — I I I 

other ellipticals: Reynolds ot al. 1 996 j- Mahadcvan 1997 ;- D4 



1988 



Matteo et al. 1999). In ADAFs, the protons and electrons 



are thermally decoupled (Coulomb collisions are rare), so 
that all the energy generated by turbulent viscous stresses 
is stored as thermal energy of the protons and ultimately 
advected beyond the black hole horizon, while the elec- 
trons remain relatively cool and radiate much less energy, 
hence the sub-Eddington luminosities. In such hot accre- 
tion flows, however, the particle mean free path is often 
comparable to the size of the system, a nd thus collectiv e 
plasma effects are likely to be significant ( |Rees et al. 1982 ) . 

In the absence of collisions, strong MHD turbulence is 
required for the angular momentum transport and energy 
dissipation in accretion flows. The way MHD turbulence 
dissipates in hot accretion flows turns out to be crucial for 
the relevance of ADAF models. Let's denote P e and P p 
the amounts of energy that heat the electrons and protons, 
respectively. In traditional ADAF models, the branching 
ratio parameter 8 = P e / P p is commonly assumed to be less 
than or similar to 1CU 2 . The predicted spectra are only 
weakly sensitive to the exact value of 5 so long as 5 <; 10~ 2 , 
while for larger values of 8 the changes in spectrum char- 
acteristics are drastic because the electrons become hot 
and radiate a large fraction of the internal energy of the 
flow via the synchrotron emission so that the flow is no 
longer advection-dominated. The above constraint on the 
value of S has been recently relaxed in the modified ADAF 



|Bisnovatyi-Kogan fc Lovelace 1997 ; |Begelman fc Chiucl: 



Narayan fc Yi 1995[ ). To proceed further, we need 



to understand what type of MHD turbulence is most likely 
to be realized in ADAFs. 

It is commonly believed that MHD turbulence in accre- 
tion flows is generated by the magneto- rotational instabil- 
ity (Balbus & Hawley 1991), which generates the large- 
scale magnetic field as well as produces random gas mo- 
tions and strong fluctuating magnetic fields in the non- 
linear regime. Hence, a variety of plasma waves may 
be generated. Only low-frequency waves may be effi- 
ciently excited by large scale motions a nd carry a signif- 
icant fraction of the turbulent energy. |Quatacrt (1998| ) 
shows that among the low- frequency MHD waves, fast 
and slow magnetosonic waves are strongly damped in the 
two-temperature T p ~ 10 12 K, T e ~ 10 9 K (T p and T e 
are the proton and electron temperatures) ADAF plasma 
by collisionlcss dissipation (known as Landau damping 
and its magnetic analog the transit time damping) while 
Alfven waves are still very weakly dam ped in the incom - 
pressible approximation. As argued by Quataert (1998), 
the "impedance mismatch" will likely inhibit excitation of 
heavily damped modes and thus the MHD turbulence in 
ADAFs will be (predominantly) Alfvenic. Additionally, it 
is natural to believe that the turbulent fluctuating mag- 
netic fields are of order of the mean, large scale magnetic 
field, B ^ Bq. and the Balbus-Hawley instability is non- 
linearly saturated. 

A detailed study particle heating processes by Alfvenic 
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turbulence in ADAFs neglecting the effects of com- 
pressibility has been perform ed by Quataert (199§| ) and 
Quataert fe Gruzinov (199S). They assumed that the 



energy which is injected into the system on large scales 
comparable to the size of the accretion flow, L ~ R, 
will cascade in fc-space to small scales, typically the pro- 
ton Larmor scale, I ~ p p = Vthp/Qp (with vth being the 
thermal particle velocity and CI — eB^/mc being the cy- 
clotron frequency) or even smaller, where it dissip ates 



Since the Alfven wave cascade is likely anisotropic (Gol 
dreicl & Sridhar 1995), the k components are related as 

fen ~ fc^ 3 L -1 / 3 . Thus, when the cascade reaches k± such 
that k±p p ~ 1, the wave energy dissipates only via Lan- 
dau damping (Cherenkov resonance) and may preferen- 
tially heat either the electrons or the protons, depending 
on the value of plasma (3 ((3 = 8ttp/Bq is the ratio of 
gas to magnetic pressure). 2 The cyclotron damping which 
may be very efficient to heat protons is unimportant in 
this case, because feii <C fcj_ and lu = kuVA <SC O p (where 
v\ = i3g/47rm p n is the Alfven speed and n is the particle 
density), so that the cyclotron resonance (u>— k\\v— fl p = 0) 
is satisfied for a negligibly small population of very ener- 
getic parti cles from the tail of the Maxw ellian distribution 
function. Quataert & Gruzinov (1999) also argue that 
only a fraction of Alfvenic energy may be dissipated on 
the scale k± ~ p p l and the rest of it passes this 'barrier'. 
On scales k± 3> p p 1 Alfven waves may be converted to 
whistler waves and cascade even further where dissipation 
on electrons may be dominant. Because the details of how 
Alfven waves are converted into whistler waves on scales 
k± ~ p^ 1 arc not known, they estimate that the transition 
from the electron dominated heating regime to the proton 
dominated one occurs somewhere inbetween ~ 5 and 
/3~ 10 2 . 

The above analysis is accurate for low-amplitude tur- 
bulence, B/B -C f, where one can neglect the effect of 
finite magnetic field pressure, B 2 /8tt, (i.e., incompressible 
plasma) and use the linear theory. In a high-amplitude 
turbulence, B / Bq ~ 1, which is likely to be present in ac- 
cretion flows, finite magnetic pressure of waves nonlinearly 
couples Alfvenic energy to ion-acoustic- like (i.e., density) 
quasi-mode perturbations, which are, in general, dissipa- 
tive. In this paper, we show that: 

1. The Alfvenic parallel cascade towards higher k\\ ex- 
ists in a compressible MHD turbulence as a result 
of nonlinear wave steepening (modulational instabil- 
ity), in addition to the perpendicular (incompress- 
ible) cascade in k±. This parallel cascade is non- 
linearly dissipative because it proceeds via the ex- 
citation of ion-acoustic (i.e., compressible) quasi- 
modes which are damped via the Landau mecha- 
nism. Even for relatively small amplitudes of tur- 
bulence, B/B ;> 0.2, and for ;> 3, much or al- 
most all magnetic energy of Alfven waves dissipates 
before it reaches the perpendicular 'dissipation bar- 
rier', k±p p ~ 1. The perpendicular cascade, thus, 
turns out to be energetically unimportant. Dissipa- 
tion in the nonlinear parallel cascade preferentially 



heats the protons and yields S <J few x 10 2 . 

2. If (3 <J 3, the turbulence bifurcates to another 'phase' 
where the Afvenic dissipation in the parallel cas- 
cade is weak and other dissipation mechanisms (at 
small scales) are required to maintain a steady- 
state. Since the parallel cascade proceeds up to 
feii ~ 4(B / Bo)(fl p /vA), the Alfvenic turbulent en- 
ergy may be efficiently dissipated due the cyclotron 
damping on protons, if ui = k\\VA — Q pi i.e., if 
B/Bq 0.3. Only protons may be heated in this 
regime. For lower amplitudes, the perpendicular 
Goldreich-Sridhar cascade is likely to dominate and 
ener gy dissipates, as discussed b y |Quataert (1998 ) 
and Quataert fc Gruzinov (1999 ). The three- wave 
cascade (Ng & Bhattacharjcc 1996), if it dominated 
the above four- wave cascade, turns out to be insignif- 
icant for electron energetics. 

We also consider other 'higher-order' nonlinear effects, 
such as particle trapping, and their relevance to particle 
heating in ADAFs. 

The paper is organized as follows. In §||we discuss the 
noisy-KNLS model of strong MHD turbulence. In §|| we 
apply the model to the conditions in ADAFs. We com- 
pare the model at hand with other competing turbulent 
processes in §[|. Section §|| is the conclusion. 

2. THE MODEL OF NONLINEAR ALFVEN WAVE 
TURBULENCE 

2.1. General considerations 

It is known that in an incompressible plasma, V • v = 0, 
the Reynolds and magnetic stresses (i.e., the fluid and 
magnetic nonlinearities) in a finite-amplitude Alfven wave 
cancel each other exactly. Such a wave behaves as if it 
is linear. Hence, only the Goldreich-Sridhar Alfvenic per- 
pendicular, k±, cascade occurs in a turbulent regime. The 
mutual cancellation of Reynolds and magnetic stresses in 
the wave breaks down when plasma is compressible. In- 
deed, the magnetic field pressure, B 2 /8ir, associated with 
the wave field exerts an additional stress onto an ionized 
gas and change its local density, n. Variations in n affect, 
in turn, the local wave phase speed, va, and intro duce a 
positive nonlinea r feedback into the wave dynamics ( pohen 
& Kulsrud 1971). This nonlinear phase speed-amplitude 



coupling in nonlinear Alfven waves results in the modu- 
lational instability which is ultimately responsible for the 
wave-front steepening, parallel, fen, Alfvenic cascade, and 
formation of collisionless shocks. The compressional non- 
linearity may also be viewed as an effective parametric cou- 
pling of a finite-amplitude Alfven wave to ion-acoustic (i.e., 
density) wave-like perturbations in a medium. 3 These ion- 
acoustic quasi-modes are compressional and have a longi- 
tudinal component of the electric field (with respect to the 
ambient magnetic field). Hence, they efficiently dissipate 
the energy via Landau damping. A quantitative theory of 
nonlinear Alfven waves is n ow well established (see a short 
review by Medvedev 199£ and references therein). Note 
that the dissipation of the finite-amplitude Alfven waves 

2 The plasma physics j3 used in this paper is different from that used in ADAF models and is related to it as /3 a< j a f = P/(J3 + 1). 
3 Such ion-acoustic modes are not plasma eigenmodes. They are driven by the Alfven wave and, hence, propagate with the Alfven speed, 
vhich is, in general, different from the sound speed in plasma. 
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is intrinsically nonlinear process a nd leads to a nu mber of 
unusual properties of such waves flMedvedcv 1999| ). 

2.2. The model 

There is only one (to our knowledge) model of strong 
nonlinear Alfven wave turbulence in a compressible plasma 
which self-consistently includes the effects of wave non- 
l inearity and particle kinetic s, e.g., Landau damping 
( Medvedev fc Diamond 1997 ). It is formulated as a 
stochastic modification of the dynamic equation of Alfven 
waves, - - the kinetic nonlinear Schrodinger equation 
(KNLS). This model of turbulence, referred to as the 
"noisy-KNLS," is thus structure-based, meaning that it 
describes the turbulence as a collection of strongly inter- 
acting coherent structures, such as shocks, solitons, cnoidal 
waves, rotational discontinuities, etc., generated by an ex- 
ternal random noise source. 4 The noisy-KNLS equation 
is 



db 



-VA 



d_ 

d.r 



(n^ 2 + N 2 UH[\b\ 2 



d 2 b 



2n p dx 2 



I (1) 



where TL is the integral operator, referred to as the Hilbert 
operator (or the Hilbert integral transform) , acting on the 
wave field: 



Here also b = (B y 



1 



|6(x')| 2 dx' 



(2) 



iB z )/B$ is the normalized complex 

wave magnetic field, \b\ 2 = \b\ 2 — (\b\ 2 ) is the fluctuating 
component of the magnetic pressure, (. . .) means appropri- 
ate averaging over the spatial domain, V means principal 
value integration, / is the random (turbulent) noise, and 
the coefficients Ni and N 2 are complicated functions of 
parameters of a plasma and are discussed below. 

The first nonlinear term in equation (nh describes the 
nonlinear (ponderomotive) coupling of finite-amplitude 
Alfven waves to ion-acoustic modes. Indeed, the first 
two terms are similar to a continuity equation for the 
wave magnetic field, d T b + d x (ub) = 0, where the "hy- 
drodynamic" velocity, u, is generated by the wave pres- 
sure, u oc tul^l 2 - The process of nonlinear steepening 
stops when it is balanced by linear Alfven wave dispersion 
which occurs at very small scales comparable to proton's 
Larmor radius, fc|| ~ p~ l , as represented by the last, lin- 
ear in b, term. The second nonlinear term in equation 
(0) describes Landau (i.e., collisionless) damping of the 
induced density perturbations, not the magnetic fluctua- 
tions. This collective kinetic effect is represented by the 
integral Hilbert operator. We emphasize that the damping 
is nonlocal which can be interpreted as the effect of par- 
ticle "memory" which is associated with particle transit 
through the wave packet. Another interesting and impor- 
tant feature is the absence of any particular scale beyond 
which the damping dominates. Indeed, the Fourier repre- 
sentation of the Hilbert operator H[f] acting on a func- 
tion / is (ife||/|fc|i \)fk = isgn(fc||)/fe, i.e., independent of 



the magnitude of . Thus, Landau damping of nonlinear 
Alfven waves occurs at all scales; it is scale invariant. We 
should comment that both nonlinearities result in the wave 
steepening and formation of sharp fronts (shocks) which 
is, practically, equivalent to the parallel cascade of wave 
energy to smaller scales. Unlike the conventional, inco- 
herent cascade paradigm, however, this cascade proceeds 
through random interactions of coherent wave structures 
(such as shocks, nonlinear waves, etc.). Hence, collective 
wave-particle interactions, such as the nonlinear Landau 
damping and particle trapping also contribute. 

We now briefly discuss basic assumptions used in the 
noisy-KNLS model. First, equation (!]) is the envelope 
equation. That is, Alfven waves themselves (carrier waves) 
are assumed to be linear and, thus, obey the dispersion 
uj = kuVA while the amplitude of these waves may vary in 
space and time and is described by equation (|l|) . The time 
variable r is the "slow" time of the large-scale dynamics of 
the wave envelope, r = (B /B) 2 tA, where tA = u^ 1 is the 
typical "Alfvenic" time. It has been assumed in deriva- 
tion that (B/Bq) 2 is small but finite. Comparison of the 
predictions of this model with observations of nonlinear 
Alfven waves in the solar wind by spacecrafts shows that 
the KNLS theory is "robust" and supports its validity even 
on the margin of applicability, (B/Bq) ~ 1. 

Second, the noisy-KNLS also assumes that waves are 
propagating in one direction, only. The existence of 
counter-propagating waves allows for another type of para- 
metric wave-wave interactions, — the decay instability, 
due to which an Alfven wave may decay into an acous- 
tic wave and a counter-propagating Alfven wave, — which 
greatly complicates any analytical treatment of the prob- 
lem. The decay instability, however, may be greatly sup- 
pressed in hot, two-temperature ADAF conditions. An 
acoustic (compressible) wave that is to be generated is 
heavily damped by collisionless dissipation. The same ar- 
gument that the "impedance mismatch" inhibits resonant 
excitation of strongly damped modes suggests here (but 
does not prove, though) that the parametric resonance dis- 
cussed above will be inefficient. 

Third, the KNLS-based theory considers planar waves, 
only. That is, no wave packet modulations are allowed in 
the plane perpendicular to the direction of wave propa- 
gation; the wave fronts are always flat. The dynamics of 
waves is, thus, purely one-dimensional. This strong limi- 
tation of the dynamics and evolution of wave structures in 
the Alfvenic turbulence has, probably, very little effect (if 
any) on the processes of particle heating. Indeed, on large 
scales, particles move along field lines, only. Thus, heat- 
ing of particles may occur (in Cherenkov resonance) via 
increase of their parallel velocity, only. The perpendicular 
wave evolution, therefore, hardly affects particle energy, 
unless k± ~ p^ 1 . Note, the effective one-dimensionality 
of the problem does not mean that waves may propagate 
along the ambient magnetic field, Bo, only. It has been 
shown that waves propagating at angle are still described 
by equation (|l]) if one formally writes the wave field as 
b = (By + B sin9 + iB z )/B , where <d is the angle 5 be- 



4 In this respect, the model is analogous to the noisy-Burgers model of hydrodynamic turbulence which describes the turbulence as a collection 
of randomly interacting collisional shocks. 

5 Strictly speaking, the angle © cannot be pushed close to 90° because dispersion becomes different for left-hand- and right-hand-polarized 
waves (dispersion on electrons vs. dispersion on protons). This sets the limit on the angle (tt/2 — 0) 2> \J m e /m p . 
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tween k and B . 

Fourth, the nonlinear coupling coefficients, Ni and N 2 , 
presented in the Appendix are calc ulated by solving the 
particle kinetic (Vlasov) equation (Spangler 1989, 1990) 
and agree wit h those calculated using the guiding cen- 
ter formalism (|Mj0lhus fc Wyller 1988|). Both methods, 



nonlinearity-to-damping coefficients : 



however, used perturbative techniques and do not include 
nonlinear effects such as particle trapping. Maxwellian 
distri bution wa s assumed for bulk particles. Gruzinov & 
Quatg !er"t (1999| ) showed that in a strong Goldreich-Sridhar 



turbulence the time-asymptotic distribution function of 
protons may be different. It is not clear whether this situa- 
tion occurs in ADAF because (i) the infall (accretion) time 
may be shorter than the velocity diffusion time required to 
establish a new distribution and (ii) the emergent distri- 
bution should, in general, be anisotropic (because the tur- 
bulence is anisotropic) and thus may be unstable against, 
for instance, fire-hose instability which will isotropize and 
thermalize the particles. Since Landau damping is a reso- 
nant process, it cares about the local slope of a distribution 
function, but its overall shape is irrelevant. Therefore, us- 
ing the Maxwellian distribution is reasonable. 

Finally, we do not specify the origin of random forcing 
for the moment. In general, involvement of a particular 
mechanism would force re-ordering in equation (|l]) and 
changes in A^i and N 2 . This complication is automat- 
ically eliminated in the renormalization-group approach, 
the key idea of which is to use bare nonlinear couplings to 
compute the noisc-rcnormalizc d (i.e., turbulent) quantities 
(Medvedev & Diamond 1997). Going ahead, the reason 
why the bifurcation point (see next section) seems to be 
determined by the bare iV's and not the turbulent ones is 
simple: both nonlinearities are cubic and renormalize in 
the same way 



2.3. The structure of Alfvenic turbulence 

We are interested in a steady-state turbulence in which 
the energy input due to random forcing is balanced by 
Landau dissipation. Since this dissipation is independent 
of \k\, it occurs at all scales; thus, no inertial interval ex- 
ists. An analytical study of the model for arbitrary noise 
is an extremely laborious task. That is why it has been 
analyzed only in the simplest case of ^- correlated in space 



and time, z ero-mean, white in k noise (Medvedev & Dia- 
mond 1997). A one-loop renormalization group technique 
has been utilized. The noise-dependent renormalized tur- 
bulent transport coefficients has been obtained. Quite im- 
portantly, they are mostly determined by properties of the 
noise source on the largest scales (i.e., the system size), — 
the so called infrared divergences. The large-fc tail of the 
noise spectrum seems to be relatively insignificant to con- 
trol global properties of turbulence. Thus, the assumption 
of white noise is likely to be unproblematic in the context 
of turbulence in hot accretion flows, because the spectrum 
of turbulence is not known, anyway. 

The existence of two distinct phases (or regimes) of com- 
pressible Alfvenic turbulence constitutes the most impor- 
tant and interesting prediction of the theory. The sys- 
tem will bifurcate (rather than smoothly transit) from one 
phase to another, as plasma parameters smoothly change. 
The bifurcation point is determined by the ratio of the 



|JVi 



\N 2 



1.3. 



(3) 



bif 



This result has been obtained rigorously from the anal- 
ysis of the self-consistency of the infrared cutoff scalings 
of the solutions for turbulent transport coefficients. The 
physical interpretation, however, is very simple and clear. 
If I A?"i I / 1 7V2 1 < 1.3, the collisionless dissipation always 
dominates the nonlinear parallel energy cascade (i.e., the 
wave steepening). Most of the energy injected on a large 
scale dissipates before it reaches the dispersion scale, and 
all the energy is dissipated via Landau damping in the 
Alfvenic cascade. Thus, the steady-state, so-called "hy- 
drodynamic" {lo, k — > 0) turbulence with no steep fronts 
is predicted. In the opposite case | JVx | / 1 JV2 j > 1.3, no 
mathematically rigorous prediction can be made, since no 
solution to renormalized equations exists. Speculatively, 
however, it is clear that the collisionless damping is weak 
compared to the nonlinear steepening. Only a small frac- 
tion of the injected energy dissipates, while most of it 
cascades along fcy to the smallest scales set by the linear 
dispersion. Since the source continuously injects energy 
into the system, it will be accumulated at the smallest 
scales. Hence, a non-steady-state, small-scale turbulence 
of Alfvenic shocklets is expected. Stationarity may be en- 
forced by including other dissipation mechanisms into the 
model. The most efficient one is the cyclotron damping on 
protons because it operates at nearly the same scales (i.e., 
the proton Larmor radius) where the wave energy is accu- 
mulated. Another mechanism which dominates for linear 
and weakly n onlinear waves (B/Bq <C 1) is the perpendic- 



ular cascade (iGoldreich & Sridhar 1995; Ng & Bhattachar 



jee 1996). The energy may cascade through k± ~ p p to 
smaller scales where Alfven waves are probably converted 
into whistl er waves. How energy is dis sipated in this case 
is unclear (Quataert & Gruzinov 1999). The discussion of 
these issues and quantitative estimates of the rates of the 
perpendicular vs. parallel cascade and cyclotron damping 
are given in 

3. APPLICATION TO HOT ACCRETION FLOWS 

3.1. Preliminary notes 

In this section, we primarily focus on the ADAF so- 
lutions which are the only known thermally stable, self- 
consistent models of hot, two-temperature accretion flows. 
The self-similar sc alings of plasma para meters re levant for 



1998): 



our p roblem are (Narayan & Yi 1995, see also Quataert 



T p ~2x 10^-JL^r- 1 K, (4a) 

T e ~ 10 9 - 10 10 K, (4b) 

(3 = const., (4c) 

B ~ lO 9 **" 1 ^! + fl-Wm-WmWr-V* G, (4d) 



Pp 
R 



g ~ 6 x 10- 9 a 1 / 2 /3- 1 / 2 m- 1 / 2 m- 1 / 2 r- 1 /4 j (4e) 



0.9(14-/3) 



-1/2--1/2 



0.37m-- 1 / 2 , 

c 



(4f) 
(4g) 
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where a is the Shakura-Sunyaev turbulent viscosity 
parameter, m — M/Mq is the mass of the central ob- 
ject in solar mass units, m = M/Ms^d is the accretion 
rate in Eddington units (AfEdd = 1-39 x 10 18 m g s" 1 ), 
and r = Rj R$ is the local radius in Schwarzschild units 
(Rs = 2.95 x 10 5 m cm). Here Bq is the strength of 
the large-scale magnetic field, (3 is assumed to be inde- 
pendent of R and takes, likely, sub-equipartition values 
(/? ;> 1 — 10), Pp/R is the proton Larmor scale in terms 
of the local scale of the accretion flow, and v r /c and 
va/c are the radial inflow velocity and the non-relativistic 
Alfven speed, 6 respectively, in units of the speed of light. 

The electron temperature is nearly constant through- 
out the accretion flow and starts to decrease at large 
radii, r J> 10 2 — 10 3 . The efficient cooling of relativis- 
tic electrons prevents higher temperatures at smaller radii. 
Thus the proton-to-electron temperature ratio ranges from 
T p /T e ~ 10 3 deeply inside the flow at R ~ few x Rs to 
T p /T e ~ 1 in the outer regions, R/R s J£ 10 3 . 

The MHD turbulence in hot, collisionless accretion flows 
is likely to originate as a result of t he magneto-rotational 
instability (Balbus & Hawlcy 1991). In the linear phase, 
the Balbus-Hawlcy instability predominantly amplifies the 
toroidal component of the large-scale magnetic field, so 
that the fluctuations are weak, B -C Bq. When the insta- 
bility reaches nonlinear amplitudes, it becomes large scale, 
ksH ~ -R -1 , so that long wavelength MHD waves are ex- 
cited and the MHD turbulence results, a s seen from numer- 
ical simulations (see, e.g., a review by Balbus & Hawley 
1998Qmd references therein). This turbulence is gener- 
ated on large scales and, hence, carries most of the grav- 
itational potential energy of the accreting gas. Since the 
ADAF gas/plasma is two-temperature, with the protons 
much hotter than the electrons, all compressional MHD 
modes are heavily damped by the collisionless (Landau 
and transit time ) damping and t hus are hard to excite, 
as suggested by puataert (1998 ). Thus, only noncom- 
pressive (Alfven) MHD modes may be efficiently excited 
and reach nonlinear amplitudes. The nonlinear satura- 
tion of the Halbus-Hawley instability means that the local 
fluctuating magnetic fields in the accretion flow become 
comparable to the averaged global magnetic field and pre- 
vent its further growth; thus in the absence of dissipation 7 
B/Bo ~ 1. At such amplitudes, noncompressibility of 
Alfven waves is no longer a good approximation and the 
nonlinear theory of strong MHD turbulence discussed in 
the previous section should be used instead. In this case, 
Landau damping will decrease the amplitude of turbulence 
to a level when dissipation balances energy input. Given 
the noise prope rties, the amplitude of turbul ence may be 
estimated from Medvedev & Diamond (1997). 

3.2. Dissipation-dominated parallel cascade 

As we discussed in the previous section (§0), there is a 
regime of strong nonlinear Alfvenic turbulence in which all 
energy that is injected into a system on large scales dissi- 
pates in the cascade and does not reach the proton Larmor 
radius scales. The bifurcation point, which is given by Eq. 

6 In the relativistic case, the displacement current cannot be neglected in Maxwell's equations. The relativistic Alfven velocity then reads 
va/^/I + vI/c 2 . We neglect this relativistic effect in our paper because it may be significant very close to the inner edge of the flow, 



(||), is a function of T p /T e and j3, only (via N± and N2, 
see Appendix). It is independent of the fluctuation level 
b oc B I Bq because steepening and damping are both cu- 
bic in b and the wave amplitude cancels out. (The rate of 
dissipation, however, does depe nd on the amplitude of tur- 
bulence, as discussed below in §|4.l].) Using the expressions 
for Ni and A2 given in Appendix and Eq. (^|), we now plot 
the borderline between the two regimes, as a function of 
plasma parameters. 

The Tp/r e -/9-diagram of state of the nonlinear Alfvenic 
turbulence is shown in Fig. |Tj. The shaded region corre- 
sponds to IA1/A2I < 1.3 and the unshaded region — to 
IA1/A2I > 1.3. The bold-faced part of the boundary curve 
represents the range of the temperature ratio relevant for 
ADAFs: T p /T e decreases from ~ 10 3 in the inner regions 
of the accretion flow at radii R <J lORs to ~ 1 in the 
outer regions at R J> 10 2 — 10 3 i?s. As is seen from the fig- 
ure, the boundary between the two phases of turbulence 
is insensitive to the temperature ratio for the ADAF con- 
ditions and set only by plasma (3. If (3 > /3 CI it — 2.6, the 
Alfvenic turbulence in an ADAF is in the regime in which 
the parallel cascade is strongly dissipative. If additionally 
the parallel cascade is dominant (i.e., faster than any other 
cascade/dissipation process), then it dissipates all injected 
energy. Alternatively, for (3 < f3 cr a — 2.6, Alfven wave en- 
ergy does not dissipate in the parallel cascade (even if it is 
dynamically dominant) and reaches k\\ ~ p~ x where it may 
be dissipated by other mechanisms, or cascade beyond this 
point to even smaller scales. 

3.3. Electron vs. proton heating in the parallel cascade 

We can also investigate the electron-to-proton heating 
ratio, 5, using the equations for Ni and N2 as follows. 
The electron and proton contributions to the wave dissipa- 
tion are proportional to t he im aginary part of the pl asma 
dispersion function, Eq. (A17). The factor in Eq. ( |A17| ) 
with integral over the distribution function has been ex- 
plicitly incorporated into the wave equation, Eq. (Q), as 
the Hilbert operator. The numerical factor was left in the 
coefficients, Eqs. ( |A6| ) - ( |A16 ), and is proportional to e~ X ' 

and e~ p for the electrons and protons, respectively. By 
vanishing by hands all those terms which are multiplied 
by e~ X " , we may easily determine how much the electrons 
contribute to the overall dissipation of waves. Calculaing 
the ratio (1 — N^/N^ = 6, where N% is the nonlinear dis- 
sipation coupling coefficient without the contribution of 
the electrons to dissipation, we obtain the heating ratio. 
Note, that in Aj the electrons still contribute to the linear 
and nonlinear dispersion of waves, as represented by the 
real part of the plasma distribution function. 

A contour plot showing the levels of constant heating 
ratio, 6 = .005, .01, .035, is superimposed onto the 
Alfvenic dominated dissipation region in Fig. |} Fig. || 
also depicts the heating ratio as a function of the tempera- 
ture ratio Tp/T e for two values of f3, f3 = 3 and 40. Clearly, 
5 depends on (3 only weakly and ranges with temperature 
from ~ .025 to ~ .005 for the typical ADAF conditions. 



»A 



near the last stable orbit, r ~ 3, only. 

7 This condition does not place any constraint on the value of Bq. 



6 



Since, the spectral ADAF models weakly depend on 6 so 
long as 6 <; few x 1CP 2 , we set a conservative upper limit 
on the electron heating in the dissipative parallel cascade 
of nonlinear Alfvenic turbulence: 



r _ Pe 

h = jr 

J- n 



£ .025, if p > p clit ~ 2.6. (5) 



4. OTHER COMPETING PROCESSES 

In the previous section (§[T^) we demonstrated that 
the energy which heats the electrons in the compressible 
Alfvenic dissipative parallel cascade in typical ADAF con- 
ditions constitues a few percents of the total energy which 
goes into the protons. If this parallel cascade would be the 
only one which operates, then equation (||) gives a solid 
prediction. In fact, there are other competing processes 
which transfer Alfvenic energy in fc-space and thus may 
affect the energetics. We consider them in order. 

In addition to the compressible parallel cascade, there 
is an incompressible perpendicular cascade due to direct 
interactions of linear Alfven wave packets. Recently two 
different theories of this process have been proposed. |Gol- 



drcich fc Sridhar (1995 ) (hereafter GS) showed that 4- 



wave interactions lead to a critically balanced anisotropic 
cascade in which fcj and k± are uniqu ely related. A 



nice paper by Ng fc Bhattacharjee (1996 ) (hereafter NB) 
demonstrates that 3-wave interactions of k\\ — pertur- 
bations are not empty and result in a more conventional 
Iroshnikov-Kraichnan-type cascade. Although there are 
some indications fro m both theoretical arguments and nu- 
merical simula tions (Goldreich & Sridhar 1997; pho & 



Vishniac 200C) that the GS cascade dominates over the 



NB cascade, the debate is still not completely resolved. For 
this reason, we consider both possibilities below. Briefly, 
we obtained that t he G S cascade puts some limitations on 
our predictions of §0. The NB cascade, however, does not 
alter our conclusions about the electron vs proton heating, 
despite that the total wave energy in the perpendicular di- 
rection may be substantial. 



4.1. Parallel 



GS 



ide 



The cri tically balanced GS cascade ( |Goldreich fc Srid- 



har 1995) is an anisotropic, (almost) perpendicular cas- 

3/2 

cade in which fcy -C k± oc k^ . This cascade involves 
linear Alfven waves and, thus, is independent of the am- 
plitude of the turbulence. No damping occurs in this pro- 
cess. The MHD (Alfven) energy is transferred from the 
scales where it is injected, k±i ~ R~ x , to the scales where 
it is dissipated or converted into other waves, kj_f ~ Pp • 
We have to compare it with the competing compressible 
parallel cascade during which much of injected energy is 
dissipated. The rate of the energy transfer depends on the 
turbulence level via the compressibility. Since these two 
cascade processes are independent (at least in the KNLS 
approximation), we may compare relative contribution of 
each of them to the electron heating as a function of the 
amplitude of MHD turbulence in ADAFs. 

The characteristic rate (i.e., the inverse e- folding time) 
of the GS cascade is the Alfvenic frequency, 7gs ~ ^>A- 
The characteristic dissipation rate in the parallel cascade 

8 Note, here we set a lower limit on fcy 
parallel cascade. 



and the nonlinear cascade rate itself are comparable and 
estimated from Eq. (pi) to be 7nl ~ ^-Vil^lH 2 . Note, al- 
though the rate of the GS cascade, 7gs oc fcy, increases 
with fcii , so does the energy dissipation rate. Thus, the to- 
tal energy dissipated in the cascade is determined by the 
width of the "inertial range," R^ 1 <; k <^ p^ 1 . The energy 
dissipation rate is, by definition, 



de 
dt 



de dk\\ 
dkn dt 



(6) 



Using the definition for the total energy cascade rate, 
dku/dt = (7gs + 7nl)^||) an d integrating the resultant 
equation, we obtain 



lne/ep = \N 2 \\b\ 2 
lnfe|| ro /fc||o |jV 2 ||6f 2 + l' 



(7) 



where eo and fc|| are initial values. The largest scale is set 
by the size of the accretion flow, feng ~ R~ 



est scale is estimated from the scaling, k» 



The small- 



of the GS cascade with kj_ m ~ P p ■ Thus, fcjj m /fc||o — 
(R/p p ) 2/3 . |Quataert fc Gruzinov (1999| ) argue that much 



of the energy dissipated at k± ~ p~ may heat electrons, 
depending on the value of (3 and numerical constants esti- 
mated from numerical simulations which are not very well 
constrained. We assume here the worst possible case in 
which all the wave energy reaching the Larmor scale heats 
electrons. 

For the electron heating to dominate by the compress- 
ible cascade, the amount of energy that reaches the Larmor 
scale must be smaller than the fraction of energy that is 
dissipated on electrons in the compressible cascade: 



:(fc||m)Ao < 5\\, 



(8) 



where 5 is given by Eq. (||). Combining all together, this 
inequality may now be re-written as follows: 



\b\ 2 > 



1 



ln<5||/ln(i?/p p ) 



\N 2 \ 1 + | In J,, / ln(i?/p p ) 



(9) 



The nonlinear coupling to dissipation, A^, depends on 
plasma parameters. The dependence of 1 7V2 1 vs. f3 for 
two values of the proton-to-electron temperature ratio: 
Tp/T e = 1, 10 3 is shown in Fig. [|. Clearly, it weakly 
depends on T p /T e and is roughly oc \f]3. To estimate the 
lower bound on |6|, we take the lowest value of | TVs | ~ 1 
for (3 ~ 3, i.e., near the bifurcation value. The normalized 
amp litude may be written in a general case as follows (see 



\b\ 2 = (B + B kJkf/B 2 Q - ((B + B Q k ± /kf/Bl) 
= {B/B a f - ({B/Bo) 2 ) + 2(B/B )k x /k 



2(B/B ) 



(10) 



for B <J Bq and an oblique and nearly perpendicular an- 
gle of propagation, k± ~ k. We take [see Eqs. @ and 

Since both cascades proceed independently, feii m may only be larger (i.e., closer to pp ) due to the 
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@)} 6 < 2.5 x 1CT 2 and p p /i? ^ 6 x 1(T 9 . Then the 
conservative limit on the fluctuation amplitude of MHD 
turbulence at which the electron heating is dominated by 
the compressible cascade is readily estimated as 



B 



£0.2. 



(11) 



Thus, for B/Bq ~ 0.2, the electron heating due to the par- 
allel, Sp and perpendicular, S±, cascades are of compara- 
ble value and the total electron heating is 6 = 6n + 5± ~ 
0.05. Hence, for the levels of turbulence which are be- 
lieved to exist in ADAFs, B/Bq ~ 1, the heating of elec- 
trons is dominated by the compressible cascade by at least 
an order of magnitude and yields 8 < 0.025, provided 
f3 > clit ~ 2.6. 



4.2. Parallel 



NB cascade 



In a similar way we consider comp etition of the NB and 
paral lel cascades. The NB cascade ( Ng Sz Bhattacharjee 



1996 ) proceeds in k± only, hence it is totally independent 
of fen . There is no dissipation during the NB cascade. The 
Alfven wave energy may be damped onl y wh en it reaches 
the k± ~ Pp 1 scale and, as discussed in §4J, will predom- 
inantly heat the electrons. To estimate the energy lost 
via the NB cascade, it is sufficient to compare the rate, 



7nb (p p 1 ) , of the NB cascade at fcj 
mum damping rate in the kn cascade, 7nl, 



to the maxi- 
ax . The latter 

is 7NL,max ~ UA, m ax\N 2 \\b\ 2 ~ v A p„ 1 \ N 2 \\b\ 2 . The former 
is estimated from Ng & Bhattacharjee (1996) as follows. 
The energy transfer rate is 



v 2 /(Nt) ~ constant, 



(12) 



where v 2 ~ E^k is the energy per unit mass, N ~ (Sv/v) 
is the number of interactions, r ~ (k^va)^ 1 is the interac- 
tion time, Sv is the perturbation, and Kyis a longitudinal 
scale of interacting Alfven wave packets. Therefore, 



7nb 



(Nr)- 1 
(E ko k )' 



e/v 2 
k 
fco 



1/2 'ft*. 

VQ 



K\\ V A , (13) 



where we used the NB scaling oc fc~ 3 / 2 and the sub- 
script "0" denotes quantities at the outer scale of the tur- 
bulence. Assuming strong turbulence, Svq ~ Vo, we obtain 
an upper limit on 7nb- Assuming also that on the outer 
scale feo ~ K o ~ K ||o ~ K ±o ~ R~ , we have 

7nb(p; 1 ) ~ tu(«o*) 1/2 x ~ «A(i?P P )" 1/2 . (14) 



Thus, the ratio is 



7nb(p p 1 ) 



7NL,: 



(WW)" (5) 



1/2 



10" 



(15) 



Thus, even if all the energy cascaded in k± heats the elec- 
trons, it constitutes about 10~ 4 of the total energy dis- 
sipated on the protons (unless B/Bq <C 1), i.e., is com- 
pletely negligible compared to equation <Sj . Note that the 
above result means only that the NB cascade is "slower" 
than the parallel cascade, while the total wave energy in 
kj_ may be non-negligible. 



4.3. Proton- cyclotron damping vs. perpendicular cascade 

The nonlinear interaction of finite-amplitude wave- 
packets re sult s in generation of high-fe harmonics, as dis- 
cussed in §2^. Such a cascade transfers wave energy along 
the local magnetic field to high fey, where the cyclotron 
resonance, to — kuv — Q p = 0, may be satisfied and the cy- 
clotron damping on (bulk) protons becomes very efficient. 
This effect was completely missing from previous studies 
because both the GS and NB cascades are (nearly) per- 
pendicular; fey is always small, so that only a very small 
number of particles from the tail of a particle distribution 
function are in the resonance. The proton-cyclotron damp- 
ing is known to heat protons only, no energy goes to elec- 
trons, since £l p <§C f2 e and electrons are off resonance. It is, 

however, difficult to rigorously estimate the \b\ 2 , because 
for the cyclotron damping to dominate, it must be faster 
than the typical GS or NB cascading time, i.e., j c J> to a- 
Such a situation may occur only when u>a — so that (i) 
Alfven waves are heavily damped and (ii) they may con- 
vert to proton cyclotron waves, which greatly complicates 
a rigorous analysis. 

One can, however, roughly estimate the fluctuation level 
at which the proton-cyclotron damping is so strong that 
most of the wave energy dissipates on protons and does 
not convert into whistlers. The parallel cascade proceeds 
until it hits a scale where the nonlinear steepening is bal- 
anced by the wave dispersion [last term in Eq. (Q)]. Thus, 
the maximum parallel wave-number to which wave energy 
cascades in the compressible cascade is readily estimated 

\b\ 2 (2n p /v A \N h2 \), where 
and 



from Eq. ([j]) to be fe|| max 

| iVi , 2 1 is the maximum of |i\Ti| and \N 2 \. The proton- 
cyclotron damping is strong when ll>a = k\\VA — ^p- Thus, 
if fc|| J> flp/vA, the MHD turbulence heats protons only. 
This occurs when 



1 



B 

~Eh ^ 4|JVi, 3 | 



(16) 



where we again took |6| 2 ~ 2(B/Bq). Here l-ZVi^l ~ 1 
for [3 ~ 3 and increases with f3 and weakly with T p /T e as 
shown in Fig. |j. For a range of (3's, 1 <; (3 <J 2.6, i.e., above 
the bifurcation threshold (dissipation in the compressible 
cascade is weak), but below equipartition, the conservative 
condition on B reads: 



B , , 

TT^O.3. (17) 



4.4. Is nonlinear trapping important? 

The parallel (coherent) cascade of the wave energy oc 
curs due to the formation and nonlinear evolution of co- 
herent wave structures, such as shocks, Alfvenic discon- 
tinuities, solitary wave-packets, etc., as discussed in §2.2 
The resonant particles, i.e., those which take the energy 
from waves and result in damping, turn out to be trapped 
in the wave potential created by the wave magnetic field. 
Those resonant particles which move a little slower than 
the wave will be reflected from the rear part of the wave 
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potential, gain energy, and start to move faster than the 
wave. The particles moving faster than the wave will be 
reflected from the front part of the wave potential and, 
correspondingly, give energy to the wave and decelerate. 
Upon a half bounce time, Tb/2, the particles which have 
been reflected from the rear of the potential cross though it 
and reach its front part, where they are reflected back and 
give their energy to the wave. The slow particles do the 
opposite at the rear part. Such a process may repeat over 
and over again, provided the wave amplitude is maintained 
constant by an external source (e.g., the Balbus-Hawley 
instability in ADAFs), to prevent particles from escape. 
Clearly, the number of slow (v < va) and fast (v > va) 
particles will oscillate with time, so that the wave damp- 
ing rate will change in time too and may become negative 
(i.e., the wave amplitude will grow) during some periods 
of time. The bounce period of a particle depends on its 
energy, because a wave trapping potential is, in general, 
anharmonic. Due to the difference in the bounce periods, 
groups of particles with different energies gain a phase shift 
which grows with time. Upon many cycles (bounces), the 
particles near the resonance, v ~ va, become completely 
randomized in phase and form an equilibrium, plateau dis- 
tribution. At this moment, the wave collisionless dissipa- 
tion quenches. 9 The process described above is referred 



we obtain 



to as the nonlinear Landau damping (see Medvedev et al 
1998[]or more details). Numerical simulations of the pro- 
cess show that it takes ~ 10 2 — 10 3 Tb or even longer for 
particles to phase mix and quench collisionless damping. 
We note here that the nonlinear Landau damping is sig- 
nificant for large amplitude waves B / Bq ~ 1, for which 
the number of trapped particles may be large. 

We now estimate whether particle trapping is important 
for the physics of ADAFs and, in particular, whether the 
Alfven wave damping during the parallel (compressible) 
cascade in ADAFs quenches. First, the kinetic energy of 
a particle in a potential is of order its potential energy, so 
that the typical velocity of a particle of mass mj (J = e,p) 
is v 3 ~ (i3 2 /87rm J n) 1 / 2 ~ VAj{B /B ). The typical scale 
of the potential, A, is set by the scale at which the energy 
is pumped into the system, i.e., the scale of the Balbus- 
Hawley instability, A ~ k 



BH 



R. Thus the characteristic 
oscillation period of a particle (the bounce time) at a given 
radial position in the accretion flow is Tb ~ A/vj, which 
yields 



T bp ~ R/(vaB/B ), T be = (m e /m p ) 1/2 r hp (18) 

for protons and electrons, respectively. Second, accreting 
gas flows into the central object and, thus, continuously 
replenished with a fresh material on the outer edge of the 
ADAF. The typical infall time at a given radius R is 



R/v r , 



(19) 



where v r is the radial velocity of a gas given by Eqs. (|J). 
Estimating the number of bounces of a particle of species j 
in an Alfven wave at a radial position R as Afj (R) ~ T T /i\,j , 



AA p ~0.4a(l + /5) 1 / 2 (S/S ), 



M e ^ 43AC 



(20) 



Note that N is independent of radius. The above equa- 
tions show that the protons in the accretion flow are always 
in the linear regime of damping (no trapping, Af p <J 1). 
The electrons may experience up to a few tens of bounces 
within the infall time, which is probably not enough to 
quench damping on electrons at all, but still may lower its 
value a little. We thus conclude that the damping process 
in the compressible parallel cascade is hardly affected by 
nonlinear particle trapping and the value of 5 ~ few x 10 -2 
represents a conservative upper limit on the the electron 
to proton heating ratio in hot, advection-dominated accre- 
tion flows with strong magnetic field turbulence, B ~ Bq. 

5. CONCLUSIONS 

The particle heating is one of the key issues for the 
advection-dominated accretion flow models. All the mod- 
els assume that protons receive most of the energy released 
due to viscous (turbulent) heating of the accretion gas, 
while electrons remain relatively cold; hence the low lu- 
minosities of ADAFs. In the collisionless plasma of the 
ADAFs, particle heating is determined solely by excita- 
tion and dissipation of plasma collective motions (waves). 
The instability of Balbus & Hawley which is believed to 
operate and produce magnetic fields in such differentially 
rotating flows may naturally result in strong MHD tur- 
bulence with highly fluctuating magnetic fields, B ~ Bq, 
during the nonlinear stages of its evolution. A detailed 
analysis of particle heating by such turbulence is presented 
in this paper. 

We show that in the most natural case of the ADAF 
parameters, (3 > 3 and B/Bq > 0.2, dissipation of the 
magnetic energy occurs predominantly in the parallel, 
compressible cascade (the effect absent for linear, low- 
amplitude Alfven waves). Most of the dissipated energy 
heats protons, while electrons receive only 6 ~ few% of 
the energy, which in agreement with the empirical values 
inferred from the spectral fits for various, low-luminosity 
accreting systems. If, alternatively, the magnetic field in 
ADAFs is close to equipartition, (3 £ 3, the energy of 
MHD turbulence is dissipated only on protons via the cy- 
clotron damping mechanism, provided the amplitude of 
the turbulence is somewhat higher, B/Bq ;> few x 10 _1 . 
The electron heating is negligible in this case. At lower 
amplitudes, however, the nonlinear effects are weak and 
the results of the linear analysis for the dissipation of the 
linear Alfven waves, addressed by other authors, hold. 
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and stimulated further work, and to Eliot Quataert for 
discussions. This work was supported by NASA grant 
NAG 5-2837 and NSF grant PHY 9507695. 



APPENDIX 

9 This process does not affect the cyclotron damping on protons, which operates at the cyclotron, and not Cherenkov, resonance. 



EXACT FORMULAE FOR THE NONLINEAR COEFFICIENTS ATj AND N 2 

The nonlinear coefficients in Eq. ([[]) has been calculated from the full kinetic treatment b y solving; the Vlasov e quation 
for arbitrary wave profile and determining the wave- induced plasma density perturbation (Spangler 1989, 199C| ). They 



are 



N 1 = M 1 -m 1 , N 2 = M 2 



m 2l 



(Al) 



where Mi and M 2 are the isotropic pressure contributions and mi and m 2 are the contributions due to pressure anisotropy 
(i.e., the temperatures and, hence, pressures may be different along the ambient magnetic field and perpendicular to it, 
P\\ ¥=P±)- 



Mi 
M 2 
mi 
m 2 



-7 I A 



AX 2 



4A2 



/e4 

fj2 



h (hh - *hh) + nf 6 (f 3 f 2 + fiU) 
ifl + nf!) 
h (hh ~ nfofi) - f 5 {hh + hh) 





(/! + 




94 (hh 


- T/3/4) 


+ Tff3 (/3/2 + /lA) 




2(/f 




93 (hh 


- TT/a/l) 


-9A(hh + hh)\ 



2 (/I + 7T/I) 



(A2) 
(A3) 
(A4) 
(A5) 



where ^/tt appears in numerators of M 2 and 7712 because the Hilbert operator carries an extra (1/ n) and additional 
(1/2) in all four coefficients absorbs it from the wave equation, as compared to Spangler (1989 ; 1990). The coefficients 
h, h, ■ ■ ■ , h, 91) • • • , 94, are defined as 



h 
h 

h 

h 

h 

h 
h 



—X p Z R (X p ) + X e Z R (X e ), 

[1 + X p Z R (X p )} + T P [1 + X e Z R (X e 

X p e- X l - X e e- X % 

X p e- X p +T r X e e- x % 
1 + X P Z R (X P ), 

X p e- X l, 
XpZn(Xp), 



(A6) 
(A7) 

(A8) 

(A9) 
(A10) 

(All) 
(A12) 



gi = [X 2 p + X z p Z R (X p ) - X p Z R (X p )} + (1/T r ) [X 2 e + X z e Z R (X e ) - X e Z R (X e )} , (A13) 

. 92 = (A p 3 - X p )er x l + (l/T r )(Xl - X e )e~ x ', (A14) 

.93 - (2A p 3 - X p )e- X l - (2Xl ~ X e )e- X * , (A15) 

.94 = [2A„ 2 + 2XlZ R (X p ) - X p Z R (X p )] - [2X 2 e + 2X 3 e Z R (X e ) - X e Z R (X e )]. (A16) 



Here Z R (X) is the real part of the plasma dispersion function for argument X 

71 



Z(X) = ±=_ r = 2* e- x * f e"^. (A17) 



The imaginary part of Z(X) appears explicitly in Eq. (|]) as the Hilbert integral operator. Physically, X = X R + iXi is 
the ratio of wave phase velocity to thermal velocity. For Alfven waves, X p and X e become: 

At last, T r = T p /T e is the proton-to-electron temperature ratio and m e /m p ~ 1/1836 is the electron-to-proton mass 
ratio. 
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Fig. 1. — A T p /T e -/3-diagram of state of nonlinear Alfvcnic turbulence. The shaded region corresponds to the compressional cascade 
dominated regime (phase) of turbulence. The contour lines represent the electron-to-proton heating ratio, <5, in MHD turbulence in ADAFs 
with B/Bo > 0.1. The unshaded region corresponds to the small-scale dissipation dominated phase, where the Goldreich-Sridhar cascade is 
dominant. The boldfaced part of the dividing (bifurcation) line corresponds to the range of the T p /T e parameter typical of ADAFs. 

Fig. 2. — The electron-to-proton heating ratio vs. T p /T e for /3 = 3 and /3 = 40. 

Fig. 3. — The coefficients of nonlinear damping, | JV2 1 (solid curves), and nonlinear steepening, N± (dashed curves), vs. j3 for two values of 
the temperature ratio, T p /T e = 10 3 (boldface curves) and T p /T e = 1 (thin curves). 
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